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Abstract. An open question of whether phenomenological fluid equations to be used 
for direct numerical simulation of turbulence are warranted on 'first principles' is 
addressed, and the problem is posed using Klimontovich microscopic density to replace 
the Boltzmann function of the classical statistical mechanics. For inert monatomic 
gases, it is shown that all the gasdynamic equations, namely, the three conservation 
equations plus the Navier-Stokes stress law and the Fourier heat conduction law 
are retrieved as governing instantaneous quantities, without having recourse to any 
concepts of averaging or statistical equilibrium. For reactive gases, however, the 
Arrhenius reaction rate law written in terms of the fluctuating temperature is not 
justified, refiecting the fact that this rate law hinges crucially on these concepts. 



1. Introduction 

Molecular approach to macroscopic behavior of gases by means of the kinetic theory 
or nonequilibrium statistical mechanics has provided firmer basis for variety of fields of 
phenomenologies. 

Most important contributions among them include; 

i) derivation of the Navier-Stokes equation from Boltzmann's equation |I|] 0, 

ii) proof that the classical 'equilibrium' thermodynamics remains exact to the first 
order in translational noneqilibrium, namely that it has the same region of validity 
as the Navier- Stokes equation, not the inviscid Euler equation and 

iii) Arrhenius' rate law derived from the Boltzmann equation for simple classical 
molecules executing reactive collisions, where the activation energy is expressed by 
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parameters of particle mechanics, namely, reduced mass and the threshold potential 
for inelastic collisions 0]. 

Hence, those widely used tools of phenomenologies are warranted on 'first principles'. 

If the flow is turbulent, however, it is not obvious that those equations remain valid 
for instantaneous flow quantities that are both stochastic and fractal In fact, as 
the matter of definition, the Boltzmann function is the mean number density in the six- 
dimensional phase(/i-)space, so is subject to statistical averaging. Therefore, the fluid 
velocity defined through a moment integral of the Boltzmann function may not resolve 
small eddies or fractal behavior. 

Therefore, if the 'first-principle' ground of the Navier-Stokes equation to be used 
for the direct numerical simulation (DNS) is asked, we have to address to unaveraged 
Boltzmann function and to its governing equation that is free from statistical concept. 
Needs for the approach along this line reflect the observation that if we leave the binary 
(long range) correlation in the Boltzmann equation, the Chapman-Enskog procedure 
applid to this equation leads not to the 'raw' Navier-Stokes equation, but to Reynolds 
averaged one 0. 

Such a function known as the microscopic density is advocated by Klimontovich in 
connection with plasma kinetic theory [|^ . It represents exact number density in the /x- 
space, a deterministic function by nature. 

The equation governing the microscopic density for a neutral ideal gas is derived in 
1^. It is formulated as an equation of continuity in the /i- space, happens to be identical 
in form with Boltzmann's equation with different implications. In what follows this 
equation will be called master Boltzmann equation (MBE). Molecular background for 
the fluid equations to be used for DNS should be a fluid version of the MBE. 

The objective of this paper is to show how these fluid equations are derived on first 
principles of molecular dynamics for monatomic inert molecules, then how the situation 
differs crucially for reactive molecules. 

2. A formalism on the microscopic density and the master Boltzmann 
equation 

Let the microscopic density be defined by the following expression after Klimontovich 


/(z,t)= E 6{z-z^^\t)) (1) 

l<s<Af 

where z = {x, v) is a point in the fi- space, namely, a 6 dimensional space consisting of 
the physical- plus molecular- velocity spaces, z^^\t) is locus of the s-th molecule in this 
space, 6 denotes the 6-D delta function and the summation spans over total number 
of molecules N under consideration. If z^^\t) is given as the solution of Hamiltonian 
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mechanics of N-body system, formula (1) provides the most accurate expression for the 
instantaneous number density in the /i- space. 

According to this 'Khmontovich formahsm' number density in the physical space is 
given by 

n{x,t) = J f{z,t)dv 

= Y: S{x-x^^\t)) (2) 
i<s<Ar 

In the same way fluid velocity u{x,t) is defined by 
nu{x,t) = j vf{z,t)dv 

= v(^\t)5ix - x^^\t)) (3) 

l<s<Af 

These expressions are to be compared with those of the Boltzmann formalism of the 
classical kinetic theory of gases; 



n 



= / fdv (4) 

Wu — j vfdv (5) 

where / is the Boltzmann function defined by 

/ = Z (6) 

In this equation overbar(~) denotes the standard 'average', namely ensumble average 
over repeated experiments, or equivalently, an average of hotwire output over a time 
long enough to smooth out fractals, yet short enough to resolve macroscopic unsteadiness 
such as periodic shedding of Karman vortices. 

We note that there are three different levels of description of the fluid velocity: 
Obviously expression (3) is the most exact one thinkable if [x^'\t) , v^'\t)] are solved 
from Hamiltonian mechanics of the N-particle system. Therefore it is a deterministic 
quantity even if the flow is turbulent. 

This quantity, however, does not correspond to the fluid velocity as observed by 
macroscopic measurements such as hotwire anemometry. Such velocity is represented 
by the one in which the microscopic density / is replaced with its spacial average 

where the integration is taken over a small, yet macroscopic volume /S.X{— Axi AX2AX3) 
such that 



nAX > 1 



(8) 
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This is where stochastic facet of turbulence makes its appearance, therefore turbulent 
fluid velocity defined in terms of / shows both stochastic and fractal behavior. 

The third level of description of fluid velocity is expression(^ that is defined by the 
Boltzmann function, where the average is taken over a long time or a large number of 
repeated experiments so that both stochastic and fractal features are averaged out. The 
third level is what is called Reynolds- average formalism or the Boltzmann formalism, 
and will not be discussed here. 

It is shown that equation for / of the first level of description has been derived as 
the equation of continuity in the /x-space. For the usual monatomic gases it is given in 
the form of Boltzmann's equation without invoking any statistical concept [§] 

'^<-D-<-§i^--l)L-m]-o (9) 

with J denoting the classical Boltzmann collision integral defined by 

Jill] = / [n^')m - mim^^^dxdv (lo) 

where the prime (') denotes a phase-space point of a particle that is supposed to reach 
the unprimed point after collision, and integration / dx spans over the physical space 
of particle z supposed to collide within unit of time. 

The master-Boltzmann equation(MBE), namely, Eq, (P) for the microscopic density, 
serves the most precise description of collective behavior of the ensemble of the N- 
particle system. It can describe thermal agitations that are not tangible to macroscopic 
sensors. See ref.pf for derivation of the Landau-Lifshitz stress/heat flux density 
fluctuation formulae with even finer temporal resolvability than the original derivation 

i- 

It should be stressed that MBE holds for /, not for / unless 

// = // 

a condition that lacks immediate justification is met. 



3. Conservation equations as fluid moments of the master Boltzmann 
equation 

Here we address to equations governing instantaneous fluid quantities by taking fluid 
moments of the MBE. A parallelism to the classical procedure by Chapman-Enskog or 
Grad [10| prevails also here except that the microscopic density replaces the Boltzmann 



function, namely that exact delta function representation replaces Hermite polynomial 
expansion around the Maxwellian distribution. 



General expression for the moment equations of the MBE (y) takes the usual form 
as follows, 



_d_ 

dXr 



VrMfdv = I{M) 



with 



I{M) = - l[[M]]f_f_dxdvdv 



M' + M' - M - M 
where M stands for one of the following 13 moment functions; 

m 2 . 1 2n 2 



m, mvj, — f;., m{wjWi SjiWf^)^ — "^k^j ("^5 mass of a molecule) (12) 



with 



Wj = Vj - Uj 



(13) 



For the first three M's that are summational invariants, I{M) is shown to vanish, and 
three conservation equations, same in form with the classical theory, result; 
dp ^ d{puj) 



dt ~^ dxi 







d d 



(14) 
(15) 

^[p(e + 1^)] + ^[piLrih + y) + + UkPj = (^6) 

Fluid quantities to appear in these equations differ from those of Chapman-Enskog 
theory in that there are no statistical concepts intervening. (See, for instance, Eqs, (2) 
and (3), in contrast with (^ and (|^).) For the rest of quantities the following definitions 
are employed; 



p^^f wlfjv = ^ E (.W - u^f5{x - x^^\t)) 

■' l<s<N 

-kT_ = — / — wifdv = -— (A:; Boltzmann's constant) 
2 Tx ./ 2 2 Tx 



(17) 



3 fc^ 

e = T 

2 m 

h = T 

2 m 



(18) 



p.^ = m j WjWifdv 



p 



m 



(19) 



l<s<7V 



6 



^ f 2 r. 



= ^ E i^'^ - ^rnvf - u,)5{x - =rW) (20) 

^ l<s<N 

where p,e,h,T,p^^ and q_. denote, respectively, pressure, specific internal energy and 
enthalpy, temperature, viscous stress tensor, and heat fiux density. All of these 
quantities correspond to those to be computed from DNS with finest (the first) level 
of resolvability, much finer than those detected by any macroscopic measurements that 
belong to the second level of description. 

4. Transport relationships for viscous stress and heat fiux density 

Conservation equations ([l^ through (|l^) are not closed; two more relationships that 



give p and q . in terms of other quantities are needed. Eq. (11) for 
— J'' — J 

1 TTl 

M- miwjWi - -Sjiwl), —wlwj 

will provide such relationships. 

For this purpose somewhat tedious calculations of 

Iji = I[m{wjWi - ^Sjiwl)] (21) 

h = ^[y^^^.] (22) 
will be carried out. Substituting (|l]) into (11) we have 

J(M) = - j [[M]]^5(a;-x(^))5(S-a;(^))(5(v-v(^))5(e-v(^~))(iSdvcie(23) 

s,s 

It is more convenient to transform velocity variables to relative and center-of-gravity 
ones as 

V = V — V, U = -{v + ■S) 1 

d{U,V) (24) 



d{v, v) 

Note also that 
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dx = J Vdn (25) 

where dfl = dQ{s, s) is the differential cross-section for a molecular encounter of a pair 
of particles (s, s) such that 

Vdn = Vbdbde = B{V, 9)dede (26) 



with b, e and 6 denoting impact parameter, azymuthal angle, and bisect angle in the 
physical space between approaching and leaving asymptote loci of s relative to s. For 
'hard' repulsive intermolecular force F = aR~^{n > 3), elementary calculation of two- 
body mechanics gives 

,2a, 2 ("-5) d3 
,V) = {-)^^VT^P-^ (27) 
m au 



where /3 is a scaled impact parameter 

HI 1 2 

p = {—)—v—b 

introduced to eliminate \^-dependence of 6'(/5); 



l-.o^-^(|)-^0 (28) 

In what follows, unless stated otherwise, we will employ Maxwellian molecules (n = 5). 
This is the particular case in which B{6,V) is independent of V, yet representing a 
reahstic model for usual neutral gas molecules. (See for example ref.[T^.) 

With this specification for the molecular property, calculation of /(M) proceeds as 
follows : 

/(M) = J BdOde 5] 5{x - x^''^)5{x - x'^'^) 

s,s 

X J ^[[MMU - ^ )5{V - w^'^ + w'^'^)dUdV 
= y 6(x - x^'^)6(x - x^'h r dt dOBie) 

xK{ ^ .w'^'^ -w^'\e,e) (29) 

where we have defined 

K{U,V,9,e)^\[[M]] (30) 
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The actual forms of K together with the integrals above for M = mwjWi and mwj^ 
have been given by Grad 0], yielding 

f -3mBi{VjVi - ^V^) 
K{U,V,9,e)B{9)d9de= I ^ 3 (31) 

[ --mB,{VjVkUk - —U,) 
respectively, where Bi is a constant defined by the following integral 

Bi = 7v [' B{9) sin^ 9 cos^ 9d9 (32) 



Eq.(31) substituted into (29) with {U, V) replaced by {\{w^'^ + w^'^), w^'^-w^''^) leads 



to 



Iji = - 3mBi ^ 6{x - x'^'^)6{x - x^'^) 



- 3mBi lY, S{x - x^'^) J2 - ( 



{s)\/ (s) (s) Sjl (s) 



+ ^ - x^'^) 6{x - x^'^){wfw\'^ - -fwf) 



(33) 



= - 2m5i 5^5(a; - x'-'^)6{x - x^'^)[^{wfwi'^' + wfwi'^') 

s,s 

( W {s? , (5) ,1/ W {s) {s) , (s) (s) (s)m 



— ABiuq 



(34) 



where definitions (0) for and (19), (20) for p.^ and q., respectively, have been made 
use of. Lines of lengthy deduction of algebra have been suppressed in deriving (34). 

The lefthand side of the moment equation (11) for the stress tensor and the heat fiux 
density are calculated using only fiuid quantities as defined here, namely, those with no 
statistical average involved, yet leading to the same expression as the 13-moment method 
(Eq. (I, 17), 0) as; 



duj dui 2 du 



and 



dxi dxj 3 dx 



5 k dT 
2 m-oxj 



(35) 



(36) 



respectively. Equating expressions (33) and (pS]) we have 



with 



^ duj ^d^_ 2^ duk. 
— dxi dxj 3 dxk' 



kT 

qbI 



(viscosity coefficient) 



(37) 
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likewise, from (34) and (p6| 



with 



b{klm) 
8Bi 



kT (heat conductivity coefficient) 



(3^ 



both exact to the ffist order smallness of and q .. 

-3 1 -3 

Formulae (37) and (38) are nothing but the Navier-Stokes and the Fourier laws 
for fluctuating fluid quantities and viscosity /heat conductivity coefficients that also 
fluctuate with temperature. 

It is rather surprising that these transport relationships that have been derived 
originally as the ffist-order deviation from the equilibrium(Maxwellian) distribution 
prove to be correct even without any statistical concept. In fact, formulae (37) and 
(38) are written in terms of collective coordinates of N-particle Hamiltonian mechanics 
alone. Their validity is subject to the only assumption that both p.^ and q. are of the 
first order smallness compared with those quantities representing ffist terms of (35) and 
(36), rspectively. 

All the quantities on the r.h.s. of (37) and (38) correspond to hotwire or 
thermocouple outputs that are stochastic and fractal. Seeming difficulty of partial 
differential representations of such fractal (therefore not differentiable) quantities are 
remedied by the experimental evidence that the fractal structure does not persist to 



infinitesimal region, but ceases to operate at Kolmogorov length |Tl|. Real difficulties 
arise when the differentials are replaced with finite differences for computational 
purposes, in which case a step size be chosen smaller than the Kolmogorov length (the 
high Reynolds number difficulty). 

5. Chemical reaction rate 

Chemical reactions of the following type 

A + B + D 

that are caused by inelastic binary collisions have reaction rates obeying Arrhenius' low 

W ~ e-^/^^ (39) 

Present has derived this rate law using a simple collision model, where the activation 
energy E is expressed in terms of molecular parameters, namely, threshold potential e 
for inelastic collisions and reduced mass m^s as 

E = ^— (40) 
ImAB ^ ' 
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It is an open question that how the rate law ( p9[ ) is modified if microscopic density 
is used instead of the Maxwelhan distributin in calculating the reaction rate; 

— " inel LBd'^AdvEVd^AB (41) 

First we note that the following relationship holds independently of the intermolecular 
force law 

cLVLab = I de I bdb 
inel. Jo Jo 

= nb^^ (42) 

where b^ denotes the maximum impact parameter for a reactive collision to occur. The 
value of b^ can be obtained from the conservation laws of angular momentum and energy 
of two-body mechanics, 



R^e = bV : const 

''^'''')iR^ + R^e') + niR) = (^)V2 - const 



where (-R(t), 0(t)) is the relative particle trajectory, Q{R) is the intermolecular potential 
as depicted in Fig.l. 

Eliminating 6 and solving for the critical impact parameter b = b'^ corr responding 
to the trajectory of the closest elastic encounter that is tangential to the ridge line 
{dR/dO = at i? = To, ^{tq) = e) of the potential 'crater', we have 

b^ , 2e ,1 



Substitution of Eqs.(42) and (g|) into (|TD gives 
W = nr' 



Vcr = (^^)^ = [2E(MX^ + M^^)]^, (M;molarweight) 



(44) 



valid for any force laws governing elastic part of collisions. Employing definition (P for 
the microscopic density, rate of reaction (44) is written in an alternative form 

s,s 

j^^^^^ 6{V - vf + v^i)V{l - ^)dV (45) 
In calculating the integral, we make use of the following formula 

/ F{V)6{V -V)dV = F{V) [ 6{V-V)dV 

JV>Vcr JV>Vcr 



FiV) J dV,dVy5{V, - %)5{Vy - Vy)S{V,; ^Vl - - V7) 
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= F{V)S{V.- ^Vl - - V^) 

= F{y)S{y-v,r) (46) 

where S is the step function defined by 

S(X;Y) = I ^' if ^ > ^ (47) 
^ ' ' 1 0, otherwise ^ ' 



Then Eq. (45) reads 



1/2 ' 



with 

y AB = -V^ 



(4^ 



This expression is the ultimate form of the instantaneous reaction rate after the 
Khmontovich formahsm; further deduction is not possible unless any assumption is 
made on Vab] its average Vab = (S/cT/vrm^^)-'^/^ for instance. Then we have 

,for kT>{l)e 
4 kT 4 (4g) 



W = a< 

where 



, otherwise 



2 / 

a = TTr^n^TiBi )• 

tttuab 



is the frequency factor. 

Another model to make formula (48) tractable would be to take the average of the 
following type; 

W = J WF{VAB)dVAB (50) 

where 

F{VAB)^-^eM-^), iclB = ^) 

is the probability density of the relative velocity distribution in statistical equilibrium. 
Then the classical Arrhenius rate law written in terms of the average temperature 

W = ae"^/^^ (51) 

is recovered as it should be. 

In Fig. 2 is shown the comparison of the two average reaction rates (49) and (pl|) 
plotted against the inverse Zeldovich number {e/kT)~^. The former which might be 
called the ignition temperature model predicts the ignition temperature Tj as 

m = (^)e 
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in somewhat overestimating manner of reality. 
6. Conclusions 

The open question of whether the fluid equations as used for the direct numerical 
simulation of turbulence, namely, those written in terms of fluctuating variables have 
'first principle' basis is looked into by the molecular approach using the microscopic 
density. On this Klimontovich formalism as contrast to the Boltzmann formalism of the 
classical statistical mechanics, three conservation equations, the Navier-Stokes viscous- 
stress law and the Fourier heat conduction law prove to be valid with much finer spatial 
resolvability, and without any concepts of averaging or statistical equilibrium. 

For reactive gases, however, the parallelism breaks down in that the Arrhenius 
reaction rate law depends crucially on these concepts, therefore that this rate law written 
by the fluctuating temperature is not justified. 
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Figure captions 

Fig.l The intermolecular potential with threshold e for reactive collisions 
Fig. 2 Arrhenius-and ignition-temperature models 



